clear;clc;
for b=1:1000
    a=sqrt((b+2)^2-b^2);
    if(a==floor(a))
        fprintf('a= %i,b= %i,c= %i\n',a,b,b+2)
    end
end

fprintf('\n')
syms u b
a=2*u+2;
b=solve((a^2+b^2-(b+2)^2),b)
c=sqrt(a^2+b^2)

b=u^2 + 2*u;
c=((u^2 + 2*u)^2 + (2*u + 2)^2)^(1/2);
for i=1:50
    a=2*i+2;
    b=i^2 + 2*i;
    c=((i^2 + 2*i)^2 + (2*i + 2)^2)^(1/2);
    if(c<=1000)
        fprintf('a= %i,b= %i,c= %i,   i= %i\n',a,b,c,i)
    end
end